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In this paper we perform a Quantal Density Functional Theory (Q-DFT) study of the 
Hydrogen molecule in its ground state. In common with traditional Kohn-Sham density 
functional theory (KS-DFT), Q-DFT transforms the interacting system as described by 
Schrodinger theory, to one of noninteracting fermions - the S system - such that the equiv- 
alent density, total energy, and ionization potential are obtained. The Q-DFT description 
of the S system is in terms of 'classical' fields and their quantal sources that are quantum- 
mechanical expectations of Hermitian operators taken with respect to the interacting and 
S system wave functions.. The sources, and hence the fields, are separately representative 
of all the many-body effects the S system must account for, viz. electron correlations due 
to the Pauli exclusion principle. Coulomb repulsion and Correlation-Kinetic effects. The 
local electron-interaction potential energy of each model fermion is the work done to move 
it in the force of a conservative effective field that is the sum of the individual fields. The 
Hartree, Pauli, Coulomb, and Correlation-Kinetic energy components of the total energy 
are also expressed in virial form in terms of the corresponding fields. The highest occupied 
eigenvalue of the S system is the negative of the ionization potential energy. The Q-DFT 
analysis of the Hydrogen molecule is performed employing the highly accurate correlated 
wave function of Kolos and Roothaan. The structure of the sources - the density, Fermi- 
Coulomb, Fermi, and Coulomb holes - as a function of the electron position are obtained, 
and from them the corresponding fields. (To our knowledge, these are the first accurate 
graphs of the Fermi-Coulomb and Coulomb holes as a function of electron position for the 
Hydrogen molecule.) As a consequence of the symmetry of the molecule, the individual fields 
- Hartree, Pauli, Coulomb, Correlation-Kinetic - are all antisymmetric about the center of 
the nuclear bond. Thus, the electron-interaction potential energy, and its Hartree, Pauli, 
Coulomb, and Correlation-Kinetic components are each symmetric about this center. The 
Coulomb correlation and Correlation-Kinetic fields, and hence their contributions to the po- 
tential and total energy are an order of magnitude smaller than those due to the Hartree and 
Pauli terms. However, the Correlation-Kinetic contribution is more significant than that due 
to Coulomb correlations. This new fact is important to the construction of approximate KS- 
DFT 'correlation' energy functionals for molecules. Finally, there is a striking similarity in 
the structure of the various sources, fields, and potential energies of the Hydrogen molecule 
for electron positions in the positive half-space encompassing one nucleus, and those of the 
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Helium atom. 
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I. INTRODUCTION 

In this paper we analyze the Hydrogen molecule {H2) in its ground-state electronic configura- 
tion (o^ls)^ from the perspective of time-independent Quantal density functional theory (Q-DFT) 
. The in principle exact framework of Q-DFT for ground and excited states, 
both nondegenerate and degenerate, has been demonstrated by application to exactly solvable 



wave functions |J,|9|. In its approximate form, Q-DFT has been applied to atoms, atomic ions, 
atoms in excited states, and positron binding, as well as to the many-electron inhomogeneity at 
metallic surfaces and metallic clusters. We refer the reader to the review articles of Refs 
for further references on these applications. This paper constitutes a first step in the application 
of Q-DFT to molecules. Here we present the essentially exact analysis of the H2 molecule via 
Q-DFT by employing the highly accurate correlated wave function of Kolos and Roothaan 11 1 . 



Beyond the understandings achieved, a principal attribute of the calculation is the knowledge 
that the structure of the corresponding Q-DFT properties for other diatomic molecules will then 
be qualitatively similar. Furthermore, these essentially exact properties can be used as the basis 
for comparison and testing of various approximations within Q-DFT prior to their application to 
more complex molecules. 

Q-DFT, in common with traditional Kohn-Sham density functional theory (KS-DFT)[l^. 
maps a system of electrons in an external field F*^^* = — Vt'(r) in their ground state to one of 
noninteracting fermions in their ground state with equivalent density. The equivalent ground-state 
energy and ionization potential are thereby also obtained. The model system of noninteracting 
fermions is referred to as the S system, S being a mnemonic for a single Slater determinant. 
(Within the framework of Q-DFT, it is also possible in principle to map into an S system in which 
the noninteracting fermions are in an excited state.) The local (multiplicative) effective potential 
energy Vs{t) of the model fermions is the sum of the external v{r) and an electron-interaction 
fee(r) potential energy, the latter being representative of all the electron correlations the S 
system must account for. These correlations are those due to the Pauli exclusion principle. 
Coulomb repulsion, and Correlation-Kinetic effects. The Correlation-Kinetic contribution is 
a consequence of the difference in the kinetic energies of the interacting and noninteracting 
systems with equivalent density. In Q-DFT, the potential energy t'ee(r) is defined as the work 
done to move a model fermion in the force of a conservative 'classical' field. The components 
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of this field each separately represent a different electron correlation. The sources of these 
component fields are quantal in that they are expectations of Hermitian operators taken with 
respect to the Schrodinger and S system wave functions. The Pauli (exchange), Coulomb 
(correlation), and Correlation-Kinetic components of the total energy are also separately expressed 
in integral virial form in terms of the fields representative of these correlations. The highest oc- 
cupied eigenvalue of the S system differential equation is the negative of the ionization potential [l^. 

The traditional KS-DFT description of the S system differs from that of Q-DFT in the 
following manner. Traditional theory is in terms of the ground state energy functional E[p\ of 
the density /o(r), and of its functional derivative. In KS-DFT, the many-body correlations noted 
above are embedded in the KS electron- interaction energy functional E^^^[p\. The corresponding 
electron-interaction potential energy of the noninteracting fermions is defined as the functional 
derivative of this functional taken at the true ground state density value: fee(r) = ^E^e^[p\/ 5p{y). 
Within KS-DFT, it is common practice to subtract the known Hartree or Coulomb self energy 
functional Eh[p\ from E(.e[p\ , thereby defining the KS 'exchange-correlation' energy functional 
E^^[p\ and its functional derivative Ua;c(r) = 5E^J^[p]/ 5p{r). The functionals {E^^^[p\,E^J^[p\) 
and their respective derivatives (uee(r), i'2;c(r)) are therefore also representative of the Pauli and 
Coulomb correlations and Correlation-Kinetic effects. KS-DFT, however does not describe how the 
different electron correlations are incorporated in the functionals {E^^^[p],E^^[p\) and hence how 
they are represented in their functional derivatives. Furthermore, the functionals {E^f,^[p\, E^i^[p\) 
are themselves unknown. As such, even if the exact wave function of an interacting system 
were known, it is not possible to construct the corresponding S system directly by following the 
prescription of KS-DFT. Hence, the potential energy Vxc{'^) is usually constructed indirectly via 
density-based methods [14-17] that employ knowledge of the 'exact' density as determined from 
ah initio calculations. 

In Section II we give a brief description of ground state Q-DFT. Section III is a description of 
the various quantal sources, fields, energies and potential energies pertaining to the S system as 
determined via Q-DFT employing the 51-parameter correlated wave function of Kolos-Roothaan. 
Concluding remarks are made in Section IV. 
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II. Q-DFT OF A NONDEGENERATE GROUND STATE 

The Schrodinger equation for a system of N electrons in an external field F^^*(r) = — Vw(r), 
and in a nondegenerate ground state, is 

[f + F + ?7]^'(X) = £;^'(X), (1) 

where T = V?, V = J2i ^(^j)' ^^^'^ ^ ~ \ Si-i |r--r.- kinetic energy, local external 

potential energy, and electron-interaction potential energy operators, ^'(X) and E are the ground 
state wave function and energy, with X = xi, X2, . . . , xjv, x = ro", and r and a the spatial and spin 
coordinates. The ground state electronic density is the expectation 

p{v) = {^\p\^) , (2) 

where p = J2i (^(r — ^i) is the Hermitian density operator. The corresponding spinless single particle 
density matrix is the expectation 

7(rr') = (^ItI^t) , (3) 

where the Hermitian operator ^ = A + iB, A = ^^j[S{rj — r)Tj{sL) + S{rj — r')rj(— a)], B = 
— |X^j[(5(rj — r)rj(a) — S{rj — r')Tj{—aL)], Tj{a.) is a translation operator, and a = r' — r. The 
diagonal matrix element of 7(rr') is the density: 7(rr) = p(r). The ground state energy is the 
expectation 

E=(^^\f+V + Lf\^)=T + Eea=t+Eee, (4) 

with the kinetic energy T = (^^\T\^^, the external potential energy Egxt = ^*|^^|*^ = 
/ p{r)v{r)dr, and the electron-interaction energy Eg^ = The ionization potential is 

/ = E^°^ — E, where E^"^ is the energy of the resulting ion when the least bound electron is 
removed to infinity. 

The differential equation for the S-system in its ground state that leads to the same density 
p(r) as that of the electrons is 

[~V^ + v{r) + veeir)](j)i{x) = eiM^); i = l,...,N, (5) 

where Vee(r) is the electron-interaction potential energy of the noninteracting fermions. The S 
system wave function is the Slater determinant ${(?!)j(x)} of the orbitals <^i(x), so that the density 
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is the expectation 

p(r) = {n<Pr}mm) = EE I'^^Mi'' (6) 

and the corresponding spinless Dirac density matrix is the expectation 

7.(rr') = {^{Mmi^i}) = J2E^ii^'')Mr'c7). (7) 

i T 



The potential energy Vge{r) is the work done to move the model fermion from a reference point 
at infinity to its position at r in the force of a conservative effective field J^^^^ (r): 

vUv) = - f :Fif{v').d\'. (8) 

The field ^^^{r) is the sum of an electron- interaction field f ee(r) representative of Pauli and 
Coulomb correlations, and a Correlation-Kinetic field Zt^{r) that is representative of the correlation 
contribution to the kinetic energy: 

:F^^/(r) = fee(r) + ^t.(r). (9) 

The field ^ee(r) is obtained via Coulomb's law from its quantal source g{rr'), the pair-correlation 
density. Thus, 

e..(r) = JS<f^,.; (10) 

where ^(rr') =< *|P(rr')|* > /p{r), and P(rr') = X^^j(5(rj — r)S{rj — r') the Hermitian pair- 
correlation operator. The pair-correlation density may be further separated into its local(static) 
and nonlocal (dynamic) components as 

girr') = p(r') + p.c(rr') (11) 
= p(r')+p.(rr') + Pc(rr'), (12) 

where the sources pxci'^r'), px{rr'), and Pc{i^r') are the Fermi-Coulomb, Fermi, and Coulomb hole 
charge distributions. The Fermi hole is defined as ^^^(rr') = —\js{rr')\'^ /2p{r), and the Coulomb 
hole is defined via Eqs.(ll) and (12). The sum rules satisfied by these charge distributions are 
/p^e(rr')dr' = -1; /p^(rr')cir' = -1; p^(rr') < 0; p^(rr) = -p(r)/2, and f Pcirr')dr' = 0. With 
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the above separation, the electron-interaction field may then be written in terms of its components 
as 

£ee{r) = eH{r) + £:cc{r) (13) 
= Sh{v)+S,{v)+Sc{v), (14) 

where the Hartree £h{t^), Pauli-Coulomb Exd'"'), Pauli Sx{r), and Coulomb £c(r) fields are due 
to their respective quantal sources p(r), pxci^r'), px{rr'), and pc{rr'). 

The Correlation-Kinetic field Zt^r) is the difference of the kinetic fields Z(r) and ^s(r) of the 
interacting and noninteracting systems, respectively: 

Zt^r) = Z,{r) - Z{r), (15) 

where Z{r) = z(r; [7])/p(r) and Zs{r) = Zs{r;[ys])/ p{r). The quantal sources of the fields 
Z{r) and Zs{r) are the single particle and Dirac density matrices. The kinetic 'force' 
z(r; [7]) is defined in terms of its components as 2;Q.(r;[7]) = 2J2/3dtcep{'r',[y])/drp, where 
*a/3(i";[7]) = \[9^/9radr'^ + d'^ /dr'i^dr'^]^{r'r")\j.i=j.i'=r is the kinetic energy tensor. The field 
Zs(r; [75]) is similarly defined in terms of the S-system tensor ta/3,s(r; [is])- 

The Hartree field Sh{t^) is conservative, and V x £h{^) = 0. This is because its source 
p{r) is a static charge, and the field may consequently be written as £h{^) = — VVFff(r), where 
Wniy) = J dr'p{r')/\r — r'\. The fields Sxcin^), ^ a;(r), and Sc{r) are in general not conservative as 
their sources are nonlocal. The sum of the fields £a;c(r) + Zt^iy) and £x{^) + ^c(r) + Zt^r) are 
always conservative. 

For systems of symmetry such that the component fields fee(r) and Zt^r) are separately 
conservative, the potential energy fee(r) may be expressed as the sum of the separate work done 
in these fields. Thus 

Vee{r) = Wee{r)+WtSr) (16) 

= WH{r) + Wxc{r) + WtM (17) 
= WH{r) + Wx{r) + Wc{r) + Wt,(r), (18) 

where Wee(r), Wnir), W^c(r), Wx{r), Wc{r), and Wt^(r) are respectively the work done in the 
fields 5ee(r), ^ if(r), £xc{r), £x{tc), ^ c(r), and Zt^{r). 
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The ground state energy is 

E = Ts + J p{r)v{r)dr + E^e + T^, (19) 

where = (^<I>{(/)j}|T'|${(/>j}^ is the kinetic energy of the noninteracting Fermions. The electron- 
interaction Eee and Correlation-Kinetic Tc energies are expressed in terms of the fields f ee(r) and 
2t^{r), respectively, in integral virial form as 

Eee = J drp(r)r • 5ee(r) and Tc = ^ J drp{r)r ■ Zt,{r). (20) 

These expressions for the energy are valid whether the fields f^ee(r) and Zt^{r) are separately 
conservative or not. Employing Eq.(13) and (14) in Eq.(20), the energy E^e may be written as a 
sum of the Hartree Eh and Pauli-Coulomb E^c (or Pauli E^ plus Coulomb Ec) energies with each 
component term expressed in integral virial form. 

Finally, the highest occupied eigenvalue of the S system differential equation is the negative of 
the ionization potential: Cm = —I- 



III. APPLICATION TO THE HYDROGEN MOLECULE 



A. Wave functions, Orbitals, and Density 



The purely electronic part of the Hamiltonian for H2 in atomic units (e = m = ?i = 1) is 

H^-lvi-lvl-±-±-^-± + l- (21) 

2 2 ria r2a rih ru 

where 1 and 2 are the electrons, and a and b are the nuclei. As the wave function of the molecule in 
its ground state is unknown, we employ the essentially exact 51-parameter correlated wave function 



of Kolos-Roothaan 



111 ] in our calculations. The symmetric spatial part of the wave function is 



vl/(rir2) = exp[-5(6 + 6)] E ^mnjkp mWivl + ^UTvfviVu (22) 

mnjkl 

with 

(i = {ria+ru)/R; ^2 = {r2a + r2b)/R; (23) 



Vi = {ria - rib)/R; r]2 = (r2a - r2b)/R, 



(24) 
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where the variational parameters are 5 and the coefficients Cmnjkpi fi2 = |ri — r2|, and 
i? = 2a is the internucleus separation. The values of the variational parameters are given in 
the Appendix. The total energy (inclusive of the internuclear potential energy Vnn = ^/R) is 
Etot{H2) = -1.174448 (a.u.) at a = 0.7005 (a.u.). The kinetic energy T = -Etot, and the total 
potential energy E^xt + -E'ee + Vun = — 2.348851(a.u.). The virial theorem ratio, which is the ratio 
of the total potential energy to twice the total energy, is 0.999981. The electron interaction energy 
component E^e = 0.58737(a.u.), and the external energy E^xt = — 3.65005(a.u.). The total energy 
[l^ of the Hydrogen molecular ion H2 at the equilibrium internuclear separation of the Hydrogen 
molecule is Etot{H2)\a=o.7005 = —0.56998 (a.u.). Thus, the ionization potential of the H2 molecule 
is / = ^tot (^2^) I a=o.7005 " ^tot(^2) = 0.60447 (a.u.). 

For two electron systems such as the Hooke's atom[l^. Helium atom, or the Hydrogen 
molecule, the orbitals of the S system in its ground (singlet) state that lead to the interacting 
system density are known. These orbitals are 4>i{v) = y^p{r)/2,i = 1,2, and are therefore known 
to the same accuracy as the wave function or density. 



The density p{0,z) along the nuclear bond z-axis is plotted in Fig.l. The density is extremely 
accurate throughout space except at and very near each nucleus. Thus, although on the scale of 
this figure, it appears that the density satisfies the electron-nucleus cusp condition[2^ exactly, in 
fact it does not. 



B. Fermi-Coulomb, Fermi, and Coulomb Holes 



For the H2 molecule in its singlet ground state, there are no correlations due to the Pauli 
exclusion principle as the two electrons have opposite spin. However, within the S system 
framework, it is customary in local effective potential energy theories to define a Fermi hole 
as px{Tr') = —p{r')/2. (This is because the pair-correlation density as determined from the 
corresponding S system wave function is g{rr') = p[r')/2.) 

In Fig. 2 we plot cross-sections through the Fermi-Coulomb />2.c(rr'), Fermi px{rr'), and 
Coulomb pc(rr') hole sources as a function of r' = (0, z') for an electron at the origin r = (0, 0) at 
the center of the nuclear bond. (Because of the cylindrical symmetry of the molecule, cylindrical 
coordinates are employed throughout.) The electron position is indicated by the arrow. The three 
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charge distributions, of course, have cyhndrical symmetry about the bond axis. More significantly, 
they are symmetrical about the electron along the z' axis. Observe that at the electron position. 



both the Fermi-Coulomb and Coulomb holes exhi 
cusp condition [2^ . (Based on the work of Ref . 



3it a cusp corresponding to the electron-electron 



2l| it is known that the wave function does not 



satisfy this cusp condition exactly. It obviously satisfies it to a good degree as evidenced by the 
figure.) As expected, at the electron position, the Fermi-Coulomb hole is more negative than the 
Fermi hole. Thus, in the region about the electron, the Coulomb hole is negative. (This is also 
the case for all the other electron positions considered.) As both the Fermi-Coulomb and Fermi 
holes satisfy the same charge conservation sum rule, there must then be regions where the former 
lies above the latter. This is clearly evident in the figure. Hence, in the outer and classically 
forbidden regions of the molecule, the Coulomb hole is positive. (The positive part of the Coulomb 
hole is more clearly evident in the figures that follow.) The Coulomb hole is both positive and 
negative as its total charge is zero. The positive part of the Coulomb hole is an indication that the 
other electron is equally likely to be in the classically forbidden region on either side of each nucleus. 



As the Fermi hole is independent of electron position, we now focus on the Fermi-Coulomb 
and Coulomb holes. In Figs. 3-5, we plot the cross-sections of these holes for electron positions at 
r = (0, a/3), r = (0, 2a/3), r = (0, a). Again, observe the cusp at the electron position for both the 
Fermi-Coulomb and Coulomb holes of each figure. Note also how the positive part of the Coulomb 
hole becomes more pronounced relative to the negative part as the electron is moved away from 
the center of the nuclear bond towards one nucleus. The positive part of the Coulomb hole is also 
largest about the other nucleus, thereby indicating that the second electron is about this nucleus. 



In Figs. 6-8, we plot the Fermi-Coulomb and Coulomb hole cross-sections for an electron in the 
classically forbidden region at r = (0,2a),r = (0,4a), and r = (0,6a). The positive part of the 
Coulomb hole continues to increase about the left nucleus at the expense of the negative part as 
the electron is moved further from the molecule. Thus, even for the asymptotic position of an 
electron at r = (0,6a), the other electron is still mainly about the left nucleus. For all electron 
positions, the Fermi-Coulomb hole />a;c(rr') is negative. 



(We note that the same cross-sections of the Fermi-Coulomb, Fermi, and Coulomb holes for an 
electron position 0.3 (a.u.) to the left of the rig ht nucleus, which corresponds approximately to our 
Fig. 4, has been plotted by Baerends et al 2^] in their study of the dissociation of the molecule. 
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However, in their figure, the electron-electron cusp in the Fermi-Coulomb and Coulomb holes is 
not present because the wave function employed by these authors is a configuration-interaction 
type wave function.) 



C. Fields, Potential Energies, and Energies. 



The electron- interaction field ifee(r), and its Hartree £h{^) and Pauli-Coulomb fxc(r) compo- 
nents along the nuclear bond axis are plotted in Fig. 9. Observe that these fields all vanish at the 
center of the bond axis or origin. This is because their corresponding sources - the pair-correlation 
density g(rr'), the density p{r), and the Fermi-Coulomb hole charge /9^c(rr') — are symmetrical 
about the center of the nuclear bond for this electron position (see Figs.l and 2). The existence 
(non-zero value) of these fields for all other electron positions is a consequence of the fact that 
their sources are not symmetrical about the electron (see Figs.l and 3-8). The fields are also all 
antisymmetric about the center of the nuclear bond. (This is a reflection of the symmetry about 
the x-y plane at the center of the nuclear bond. As such the potential energies obtained from these 
fields will be symmetric about this point.) In the positive half-space, there is a maximum in the 
electron-interaction and Hartree fields, and a minimum in the Pauli-Coulomb field. The Hartree 
and Pauli-Coulomb fields are of the same order of magnitude and opposite in sign. This is because 
their sources, p(r) and /9a;c(rr') respectively, are of the same order of magnitude and opposite in 
sign. Asymptotically, in the z direction these fields decay as Seei^) ~ l/-^^, £h{y) ~ and 
£'xc(r) ~ as they must . (It is interesting to note that with a slight translation to the 



right, theplots of the fields in the positive half-space, are strikingly similar to those of the Helium 



atom 



The Pauli ^'^^(r) and Coulomb £'c(r) field components of the Pauli-Coulomb field £'xc(r) along 
the nuclear bond axis are plotted in Fig. 10. Again, these fields vanish at the origin and are 
antisymmetric about it. Hence, the corresponding potential energies obtained from these fields 
will be symmetric. In the positive half space, the Pauli field £'x(r) is negative as its source is 
a negative charge. The Coulomb field £'c(r), on the other hand, is positive in the inter-nuclear 
region and negative throughout the region beyond the right nucleus. This structure is attributable 
to the fact that the Coulomb hole has both a positive and negative component. Asymptotically, 
the Pauli field decays as fa;(r) —l/z^, whereas the Coulomb field £'c(r) has essentially vanished 
by about z = b (a.u.). (Once again in the positive half-space, the structure of these fields when 
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translated slightly to the right, is similar to those of the Helium atom. In particular, we note that 
the structure of the Coulomb holes of the Hydrogen molecule for electron positions z > a (see 
Figs. 6-8) is vervsimilar to those of the Helium atom for electron positions away from its nucleus 
(see Figs. 3,4 of |9|).) As is the case for atoms, it turns out that the asymptotic structure along 
the nuclear bond axis of {jr^-'^-^(r) — £'/f(r)} ~ <S'x(r) ~ — 1/z^ . Thus, the asymptotic structure 
of the electron-interaction potential energy t'ee(r) minus the Hartree potential energy Wh{t^) is 
again due to Pauli correlations: {t>ee(r) — Wnir)} ^ Wx{r) ~ —^/z as shown in Fig. 11. 



Since in the S system description of two electron systems fa;(r) = — i5'//(r)/2, the curl of the 
Fermi field along the nuclear bond z axis direction vanishes: V x ^'^.(r)!^ = 0, as it does in all 
directions. Hence, the work done Wx{0,z) plotted in Fig. 11 is path independent. Along the 
nuclear bond axis, however, the V x £ci^)\z ^ and V x Zt^{r)\z 7^ . But in this and all 
directions, the curl of the sum of the fields fc(r) and Zt^{r) vanishes: V x [6'c(r) + ^ic(r)])|2 = 0. 
Therefore, the work done in the sum of these fields in all directions, and hence along the nuclear 
bond axis Vc{0,z) = Wc{0,z) + Wt^{0,z) is path independent. The calculation of the potential 
energy fc(0, z) is straightforward. However, our use of the Kolos-Roothaan wave function, in spite 
of its accuracy, leads to Vc{0,z) being singular at the nucleus. This occurs due to the component 
Wt^ (0, z) that requires a cancellation of the kinetic fields of the interacting and noninteracting 
systems. The underlying reason for the singularity, however, is that the wave function does not 



satisfy the electron-nucleus cusp condition exactly. In a recent paper [23|, we have proved by 
employing the integral form of the electron- nucleus cusp condition that in local effective 
potential energy theories and for arbitrary symmetry, the potential energy i'ee(r) is finite at the 
nucleus. Furthermore, it is shown that this finiteness is a direct consequence of the satisfaction 
of the electron-nucleus cusp condition by the Schrodinger wave function. (As a consequence, for 
example, this potential energy is singular at each nucleus when determined either from Gaussian 
geminal [2^ or configuration interaction 25] wave functions.) Hence, in order to obtain Vc{0,z), 
we have employed our calculated results in regions other than near the nucleus, and smoothed 
the curve through each nucleus. ( A comparison of our results with the work of Gritsenko et al 3] 
who in their self-consistent calculations assumed Wee(r) to be finite at the nucleus show the two 
curves to be indistinguishable throughout space.) The potential energy Vc{Q,z) is plotted in 
Fig. 12. Observe that fc(0, z), and thus the sum of the Coulomb and Correlation-Kinetic potential 
energies is an order of magnitude smaller than Wx{Q,z), the Pauli contribution. The potential 
energy Vc{^,z) has considerable structure, is symmetric about the origin, and is mainly positive. 
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indicating thereby that its principal contribution is Correlation-Kinetic. (Recall that the Coulomb 
field is principally negative in the right-half space (see Fig. 10) so that the Coulomb potential 
energy TVc(r) is negative.) The plot of Vc{0,z) translated to the right nucleus is very similar in 
shape and magnitude to the corresponding potential energy fc(r) of the Helium atom (see Fig. 4 

of H). 



To obtain a quantitative sense of the separate Coulomb and Correlation-Kinetic contributions 
to Vc{0,z), we plot in both Figs. 11 and 12 the work done Wc{0,z) along the path of the nuclear 
bond in the force of the Coulomb field £^c(r). From Vc{0,z) and Wc{0,z) we obtain Wt,.{0,z) 
which is also plotted in Fig. 12. The corresponding Correlation-Kinetic field Zt^{0,z) is shown 
in Fig. 13. Note that this field too is antisymmetric about the origin. Once again, there is a 
striking similarity between the plots of Wc{0,z), Zt^{0,z), and Wt^{0,z) when translated to the 
right nucleus to those of the corresponding properties of the Helium atom The Coulomb 
correlation part WciO,z) is negative throughout space and vanishes by about z = 5 (a.u.). 
The Correlation-Kinetic piece Wt^ (0, z) is throughout positive and asymptotically decays more 
slowly. The field Zt^{0,z) is principally positive throughout space. Thus, the Correlation-Kinetic 
energy is positive: T = 1.1745 (a.u.), = 1.1414 (a.u.), = 0.0331 (a.u.) (The cor- 
responding value of Tc for the Helium atom is 0.0365 (a.u.) 9]). (We note that the Wc{0,z) 
and Wt^(0, z) do not each separately represent a potential energy. Their sum which is Vc{0, z) does.) 



IV. CONCLUDING REMARKS 



This paper is the first application of the Q-DFT quantal source and field perspective to a 
molecule, and much has been learned as explained in the previous section. The symmetry of 
the H2 molecule dictates that the individual fields £'x(r), £h{^)-, ^c(r), 2t^{r) representative of 
the Pauli and Coulomb correlations, and Correlation-Kinetic effects respectively, must each be 
antisymmetric about the center of the nuclear bond. The corresponding electron-interaction 
potential energy Veei^) representative of these correlations as determined by the work done in 
the force of these fields is then symmetric about this point as also dictated by the symmetry of 
the molecule. The potential energy Vee(r) is also finite at each nucleus, as must be the case |23l |. 
The Hartree ^//(r) and Pauli ^'^.(r) fields are the largest in magnitude and opposite in sign, the 
former being positive and twice as large as the latter. As such the principal contributions to 
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the electron-interaction energy Eee and potential energy Vee(r) are due to the Hartree and the 
Pauli correlation terms. The Coulomb fc(r) and Correlation-Kinetic 2t^{r) fields tend to cancel 
each other, so that the contribution of their sum to the potential energy Wee(r) is an order of 
magnitude smaller. However, as the potential energy component fc(r) representing the sum of 
these correlations is principally positive (see Fig. 12), it is evident that the Correlation-Kinetic 
effects are more significant. They are also more significant asymptotically, where the Coulomb 
correlation contributions to the potential energy vanish. Thus, Correlation-Kinetic effects play an 
important role in local effective potential energy theories of the H2 molecule. We further note 
that in the construction of approximate KS-DFT 'exchange-correlation' and correlation energy 
functionals i?^'^[r] and E^'^[r] for molecules, Correlation-Kinetic effects must be incorporated if 
an accurate S system representation of molecules is to be obtained. 



On the basis of the Q-DFT results determined from the H2 molecule, it is evident that 
the qualitative features of the quantal sources, fields, and potential energies for other diatomic 
molecules will be similar. However, the fields and hence the potential energies of these diatomics 
will have more structure as a consequence of the additional molecular subshells. We expect this 
added structure to be similar to that observed in atoms as the number of shells is increased. 
Finally, we note that the accuracy of approximation methods within Q-DFT and KS-DFT 
can be tested by comparison with these essentially exact results. 



We conclude by reiterating the striking similarity between the Q-DFT properties of the 
Hydrogen molecule and the Helium atom for electron positions in the positive half space. It is 
interesting that in spite of the presence of a second nucleus, and therefore of a different symmetry, 
the quantal sources and fields representative of the various electron correlations in the Hydrogen 
molecule are so similar to those of the Helium atom. This speaks to the commonality of properties 
of these distinct quantum systems as exhibited within the framework of Q-DFT. 
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TABLE I: Variational parameters in the normalized 51-parameter correlated wave function for the ground 
state of i/alli- 
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APPENDIX: Wave function parameters 

The values of the parameter b and the coefficients Cmnjkp for the wave function of Eq.(21) are listed 
in the table I. 
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0.5 




z{a.u.) 

FIG. 1: The electron density p(0, z) of the hydrogen molecule along the nuclear bond axis in atomic units 
(a.u.). The nuclei are on the axis at a = 0.7005 (a.u.) indicated by the two dots. 
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FIG. 2: Cross-sections of the Fermi-Coulomb pa;c(rr'), Fermi px{rr'), and Coulomb Pc(rr') holes along the 
nuclear bond axis for an electron at the center r = (0, 0) of the bond. The electron position is indicated by 
the arrow. 
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FIG. 3: Cross-sections of the Fermi-Coulomb /9xc(rr') and Coulomb Pc(rr') holes along the nuclear bond 
axis for an electron at the center r = (0, a/3) of the bond with the electron position indicated by the arrow. 
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FIG. 4: The same as in Fig. 3, but with the electron at r = (0, 2a/3). 




FIG. 5: The same as in Fig. 3, but with the electron at r = (0, a). 
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FIG. 6: The same as in Fig. 3, but with the electron at r = (0, 2a). 
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FIG. 7: The same as in Fig. 3, but with the electron at r = (0,4a). 
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FIG. 8: The same as in Fig. 3, but with the electron at r = (0, 6a). 
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FIG. 9: The electron-interaction £ee{0,z) field, and its Hartree £h{0,z) and Pauli-Coulomb £xc{0,z) 
components along the nuclear bond axis. 
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FIG. 10: The Pauli £x{0,z) and Coulomb £c{0,z) fields along the nuclear bond axis. The function —l/z^ 
is also plotted. 
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FIG. 11: The Pauli potential energy Wx(0, z) along the nuclear bond axis. The work done Wc(0, z) in this 
direction in the force of the Coulomb field £c(0, z), and the function —1/z, are also plotted. 
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FIG. 12: The potential energy Vc{0, z), the sum of the Coulomb and Correlation-Kinetic potential energies, 
along the nuclear bond axis. The work done Wc(0, z) of Fig. 11, and the work done Wt^(0, z) in the force 
of the Correlation-Kinetic field Zt^{0, z) of Fig. 13, are also plotted. 
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FIG. 13: The Correlation-Kinetic field (0, z) along the nuclear bond axis. 



